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ABSTRACT 


This  report  summarizes  the  work  performed  under  the  contract 
N60530-79-C-0006  for  the  period  from  October  23,  1978  to  September 
30,  1979.  It  Includes  the  study  of  dynamic  burning  effects  In  the 
combustion  of  solid  propellants  with  cracks,  and  the  use  of  granular 
bed  combustion  codes  for  solving  hazard  problems  In  rocket  propulsion 
systems. 

The  Zeldovich  quasi-steady  flame  model  was  used  In  the  evaluation 
of  the  dynamic  burning  effect.  The  dynamic  burning  rate  was  obtained 
by  solving  the  transient  heat  conduction  equation  for  the  solid 
propellant  and  using  the  Zeldovich  map  for  determining  the  heat 
feedback.  Results  indicate  a  stronger  dynamic  burning  response  for 
a  higher  pressurization  rate,  larger  energy  storage  in  the  propellant, 
and  lower  Initial  pressure.  A  dynamic  burning-rate  augmentation  function 
was  developed  to  facilitate  the  Incorporation  of  the  transient  burning 
effect  into  the  crack  combustion  code.  The  augmentation  function  is 
In  close  correlation  with  results  obtained  from  the  finite  difference 
method  over  the  broad  range  of  conditions  studied. 

The  mobile  and  fixed  granular  bed  combustion  codes  (MGBC  and  FGBC), 
users  manuals  for  these  two  programs,  sample  input  data,  and  output 
listing  were  delivered  to  NWC. 
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NOMENCLATURE 


SYMBOLS 


*1»  a2’  a3»  at 
A 
BR 


bl‘  V  b3 

cl*  °2 *  C3 


c 

C 

E 


Dimensionless  coefficients  used  in  Eq.  (32) 

Coefficient  of  the  Arrhenius  Expression.  Eq.  (7) 

Burning  rate  augmentation  parameter  defined  in  Eq.  (30) 

Coefficients  defined  in  Eq.  (13) 

Coefficients  defined  in  Eq.  (14) 

Specific  heat  at  constant  pressure 

Dimensionless  parameter  used  in  Eq.  (24) 

Activation  energy  for  surface  reactions 

A  dimensionless  factor  used  to  Increase  the  mesh  size  in 
the  y  direction 

Coefficients  used  in  Eq.  (33) 

Pressure 

Initial  pressure 

Dimensionless  initial  pressure  defined  in  Eq.  (31) 
Dimensionless  pressurization  rate  defined  in  Eq.  (28) 
Burning  rate  cf  the  propellant 
Universal  gas  constant 
Time 

Temperature 

Vertical  distance  measured  from  the  propellant  surface 
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NOMENCLATURE  Cont'd. 


SYMBOLS 


a 

X 

P 


T 

P 

T 

8 


Thermal  diffusivity 
Thermal  Conductivity 
Density 

Propellant  temperature  sensitivity 

Dimensionless  time  defined  in  Eq.  (29) 

Characteristic  time  of  the  gaseous  flame 

Transient  pressure  variation  time  (characteristic 
time  of  the  pressure  variation) 

Characteristic  time  of  the  unburned  solid  phase 

The  temperature  gradient  at  the  propellant  surface 


SUBSCRIPTS 

g  Gas 

i  i1*1  nodal  point  as  shewn  in  Fig.  1,  or  initial  condition 

pr  Propellant 

ps  Propellant  surface 

SUPERSCRIPTS 

j  jth  time  step 

*  reference  quantities 

o  steady-state  conditions 
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1.  INTRODUCTION 

1 . 1  Dynamic  Burning  E(f«ct»  In  the  Combustion  of  Solid  Propellant  Cracks 

1.1.1  Background 

In  solid  rocket  propulsion  systems,  the  Instantaneous  burning 
rate  of  a  solid  propellant  under  rapidly  changing  non-steady  pressure 
conditions  differs  greatly  from  the  steady-state  burning  rate  value 
corresponding  to  the  same  pressure.  In  general,  dynamic  burning  13  most 
frequently  encountered  under  highly  transient  conditions,  especially 
during  the  ignition  and  extinction  phases  of  combustion  processes.  The 
deviant  burning  behavior  that  occurs  during  these  periods  alters  the 
combustion  chamber  dynamics  and  flame-spreading  rates,  thus  directly 
Influencing  system  performance.  For  example,  both  the  thrust  level  and 
the  burning  time  of  a  high-performance  solid  rocket  may  be  significantly 
altered  because  they  depend  upon  the  burning  rate  of  the  propellant. 

The  dynamic  burning  behavior  of  solid  propellants  has  been 
Investigated  by  various  researchers,  both  theoretically  and  experimentally. 
An  extensive  literature  review  of  the  research  work  on  dynamic  burning  Is 
given  In  Ref.  1. 

Physically,  the  dynamic  burning  effect  is  introduced  during 
rapid  pressure  changes  caused  by  the  finite  time  interval  required  for 
the  temperature  profile  of  the  condensed  phase  to  follow  transient 
pressure  variations  (Ref.  2).  Under  an  extremely  rapid  pressure  excursion, 
the  temperature  profiles  In  the  gaseous-flame  and  surface-reaction  zones 
may  also  lag  behind  the  pressure  variation.  The  deviation  of  the  transient 
burning  rate  from  the  steady-state  burning  rate  depends  upon  two  opposing 
phenomena:  (a)  the  preheating  effect  of  the  propellant  surface,  and  (b) 

the  out-of-phase  blowing  effect  of  the  chemically  reacting  gases  adjacent 
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to  the  burning  surface.  The  net  result  of  these  phenomena  nay  result  in  a 
burning-rate  overshoot.  Experiments  indicate  that,  in  addition  to  these 
tiro  primary  effects,  there  are  often  numerous  secondary  effects  which  may 
Influence  the  heat  feedback  to  the  propellant  surface  and  therefore  alter 
the  instantaneous  burning  rate*  Secondary  effects  may  include  change  of 
chemical  kinetics  in  the  flame  zone,  change  of  flame  luminosity,  change  of 
the  fraction  of  Instantaneous  heat  fluxes  to  fuel  binder  and  oxldltar 
surfaces,  etc.  These  interrelated  physical  phenomena  are  illustrated  in 
the  flow  chart  given  in  Ref.  1. 

The  conditions  under  which  the  effect  of  dynamic  burning  is 
Important  are:  (1)  if  the  characteristic  times  associated  with  the 
gaseous  flame,  T  ,  and  the  solid  propellant,  T  ,  are  not  negligibly  small 
in  comparison  with  that  of  the  transient  pressure  variation  time,  T  ;  or 
(2)  if  the  thermal  wave  penetration  depth  is  large  enough  for  the  Initial 
preheating  effects  to  introduce  a  significant  burning-rate  augmentation 
effect.  In  the  convective  burning  of  a  solid  propellant  with  initially 
open  or  submerged  cracks,  the  pressure  variation  time  can  be  extremely 
short,  especially  during  the  transition  from  convection  burning  to  detona¬ 
tion.  Therefore,  the  dynamic  burning  effect  could  be  important  in  the 
flame-spreading  and  crack-combustion  processes. 

1.1.2  Motivation  and  Objectives 

In  the  theoretical  model  developed  at  The  Pennsylvania  State 
3-5 

University  to  describe  the  combustion  processes  in  solid-propellant 
cracks,  the  effect  of  dynamic  burning  has  not  yet  been  taken  into  account. 
However,  under  certain  conditions  the  dynamic  burning  effects  cannot  be 
Ignored.  More  specifically,  the  dynamic  burning  behavior  becomes  important 
when  the  characteristic  time  associated  with  the  unburnt  condensed  phase  is 
of  the  same  order  of  magnitude  or  greater  than  that  associated  with  the 
transient  pressure  variation.  The  investigation  of  the  dynamic  burning 
behavior  in  a  solid  propellant  crack  is,  therefore,  of  considerable 
Importance. 
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Tne  objective  of  the  dynamic  burning  study  are  as  follows: 

1.  To  Improve  the  crack  combustion  code  to  account  for  the 
dynamic  burning  by  incorporating  a  suitable  dynamic  burning 
model. 

2.  To  evaluate  the  significance  of  dynamic  burning  by 
conducting  *  parametric  study  for  different  pressurization 
rates,  initial  pressure,  and  initial  subsurface  temperature 
profiles. 

3.  To  develop  a  burning-rate  augmentation  function  for 
efficient  computation,  and  to  include  it  into  the  crack 
combustion  program. 

1 . 2  Literature  Survey  of  Recent  Studies  in  Granular  Propellant  Bed 
Combustion 


Combustion  of  granular  propellants  jnder  a  strong  confinement  may 
lead  to  the  transition  from  deflagration  to  detonation.  A  number  of 
models  have  been  proposed  to  simulate  the  processes  of  granular  propellant 
combustion.  These  models  were  developed  for  the  study  of  gun  interior 
ballistics.  A  comparison  of  the  theoretical  models,  and  experimental 
measurement  techniques  is  summarized  in  Ref.  6. 

It  lr  the  purpose  of  the  literature  survey  to  summarize  additional 

research  investigations  madr.  in  this  area,  following  the  1975  JANNAF  work- 
6 

shop.  The  literature  surxey  is  aimed  particularly  at  the  applicability 
of  existing  theoretical  models  for  the  simulation  of  rocket  propulsion 
problems  associated  with  the  deflagration-to-deionation  transition  (DDT). 
The  survey  will  also  aid  in  determining  the  technological  gaps  and 
Inadequacies  of  existing  models  used  in  DDT  studies. 
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1.3  Delivery  of  Fixed  and  Mobile  Granular  Bed  Combustion  Codes  (FGBC  and 

MGBC)  to  NWC  and  Consultation  for  the  Implementation  of  the  Two 

Prog raws 

In  recent  years  the  Principal  Investigator  has  developed  two  granular 
bed  combustion  (GBC)  models.  One  was  designed  for  a  fixed  bed  and  the 
other  for  a  mobile  bed.  The  mobile  granular  bed  combustion  (MGBC)  model 
is  more  general  than  the  fixed  granular  bed  combustion  (FGBC)  model  since 
the  motion  of  the  propellant  grain  is  considered  In  MGBC.  However,  the 
computation  time  for  the  mobile  bed  Is  considerably  longer  in  comparison 
to  that  of  the  fixed  bed.  Depending  upon  the  nature  of  the  problem,  one 
model  may  have  advantages  over  the  other.  These  programs  have  been 
successfully  executed  on  the  IBM  370  computer,  but  may  require  modifica¬ 
tions  for  other  types  of  computers.  The  users  at  NWC  may  also  need  a 
certain  amount  of  help  and  consultation  In  implementing  these  programs. 
Since  written  documentation  of  these  programs  was  not  available,  manuals 
for  users  of  the  programs  should  be  prepared. 
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H.  DESCRIPTION  OF  WORK  PERFORMED 


2. 1  Dynamic  Burning  Effects  In  Combustion  of  Solid  Propellant  Cracks 

2.1.1  Theoretical  Model 

As  pointed  out  In  the  Introduction,  the  dynamic  burning  model 
takes  Into  consideration  the  time  lag  for  the  variation  of  the  temperature 
profile  In  the  solid  propellant  under  a  pressure  excursion.  This  time  lag 
manifests  Itself  In  a  time-dependent  burning,  called  the  dynamic  burning.  In 
general,  the  dynamic  burning  rate  is  a  function  of  the  pressurization  rate, 
local  pressure,  and  the  physiocheraical  properties  of  the  propellant.  The 
following  assumptions  have  been  made  in  the  analysis  for  mathematical 
tractablllty  of  the  problem: 

1.  The  gas-phase  reaction  zone  is  considered  to  be  quasi-steady 
In  other  words,  because  the  relaxation  time  associated 

with  the  flame  is  much  shorter  than  that  associated  with 
transient  pressure  variation,  the  flame  adjusts  itself 
Immediately  to  chamber  conditions. 

2.  No  subsurface  chemical  reactions  are  present.  All  chemical 
reactions  are  confined  to  the  propellant  surface  or  in  the 
gas  phase. 

3.  Propellant  surface  temperature  is  uniform  and  the 
propellant  Is  homogeneous. 

4.  Cross-flow  velocities  at  the  propellant  surface  are  small; 
therefore,  no  simultaneous  erosive-burning  effect  is  present. 

The  validity  of . assumption  1  can  be  shown  by  comparing  the 
characteristic  times  associated  with  the  pressure  variation,  the  condensed 
phase,  and  the  gaseous  flame  zone.  In  the  crack  combustion  experiments,  the 
typical  time  associated  with  the  pressure  excursion  is  in  the  order  of 
1  ms.  The  two  other  characteristic  times  are  evaluated  as  follows: 
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2 

t  -  a  / r.  3  1.8  ms  for  r,  ■  1  cm/s 

s  pr  b  b 


(1) 


T  «  —  T  a  0.01  T  ■  0.018  ms 

8  Wls  s 


(2) 


Since  <<  r the  quasi-steady  flame  assumption  is  valid.  It  Is  also 
clear  that  Tg  and  are  of  the  same  order  of  magnitude;  therefore,  the 
dynamic  burning  effect  must  be  considered. 


Assumption  2  Is  reasonable  fcr  AP-based  composite  solid  pro¬ 
pellants  In  the  pressure  range  of  Interest.  Even  though  some  propellants 
may  have  subsurface  heat  release,  the  heat  release  zone  Is  so  thin  and  so 
close  to  the  surface  that  any  subsurface  heat  release  can  be  lumped  with 
surface  heat  release. 


Assumption  3  Is  employed  for  mathematical  simplicity.  It 
allows  the  use  of  a  one-dimensional  transient  heat  conduction  equation 
for  the  solid.  Researchers  In  the  field  have  not  as  yet  Investigated  the 
three-dimensional  treatment  of  heat  conduction  processes  In  composite 
solid  propellants  because  of  the  complexities  associated  with  the  three- 
dimensional  structure  of  oxidizer  crystals  and  fuel  binder  and  the  fact 
that  numerical  solutions  are  cumbersome  and  time  consuming. 

Assumption  A,  that  there  is  no  erosive  burning  effect,  is  not 
valid  near  the  crack  opening  region.  However,  this  assumption  does  not 
necessarily  introduce  significant  error  in  the  overall  prediction  of  the 
combustion  processes  in  cracks,  especially  for  long  cracks  with  high 
degree  of  tortuosity. 

Methods  generally  used  for  solving  the  quasi-steady  flame 
model  can  be  classified  into  the  following  categories: 
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1.  dP/dt  method 

2.  Flame  description  metl  od 

3.  Zeldovich  method 

The  dP/dt  approach  is  restricted  to  those  cases  which  deviate  only 
slightly  from  the  steady  state,  and  thus  is  not  suitable  for  the  crack 
combustion  problem.  The  flame  description  approach  is  inconvenient 
since  it  often  requires  detailed  knowledge  of  the  flame  structure  and 
propellant  surface  reaction  phenomena. * 

The  Zeldovich  approach  appears  to  be  most  suitable  for  the 
study  of  dynamic  burning  in  solid  propellant  cracks.  In  this  model,  the 
uncertainties  of  modelling  a  solid  propellant  flame  are  bypassed,  and 
the  measurable  steady-state  burning  characteristics  are  used  to  deduce 

the  non-steady  heat  feedback.  This  method  can  be  used  as  long  as  the 

■) 

quasi-steady  flame  assumption  is  valid.” 


The  governing  equation  for  the  solid  propellant  is  provided 
by  the  following  one-dimensional  transient  heat  conduction  equation: 


Pr  -  r  _Pl  +  a  _ El 

3t  b  ay  pr  ay2  (3) 

where  y  is  measured  from  the  instantaneous  propellant  crack  surface  at  any 
axial  location  as  shown  in  Fig.  1. 


The  initial  and  boundary  conditions  are: 
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FIG. 


Flow  direction  along 
the  crack  channel  - - ► 


T 

pr 


(1) 


T 

ps 


Description  of  the  physical  model. 
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Before  Ignition  of  the  propellant  at  a  given  axial  location,  the  burning 
rate  r^,  in  Eq.  (3)  is  zero  and  the  subsurface  temperature  gradient 
Is  obtained  from  a  well-established  empirical  correlation  (Dittus-Boelter* s 
correlation)  for  convective  heat  transfer  In  a  flow  channel. 

After  the  onset  of  ignition,  the  Instantaneous  burning  rate 
Is  related  to  the  propellant  surface  temperature  by  an  Arrhenius  pyrolysis 
expression: 


r.(0  ■  A  e  RuTps  (7) 

b 

The  subsurface  temperature  gradient  In  a  functional  form  can  be  written 
as 


♦g(t)  “  ♦8trb(t>» 


(8) 


Following  Zeldovlch's  concept  of  quasi-steady  flame,  this 
function  Is  the  same  as  that  obtained  from  steady-state  burning  conditions. 
A  schematic  drawing  showing  the  method  for  obtaining  the  burning  rate  map 
[♦  ■  ♦  (r^,  from  the  steady-state  experimental  data  is  shown  In  Fig. 

2.  Essentially,  the  dynamic  burning-rate  map  Is  constructed  through  the 
use  of  (a)  the  pyrolysis  data,  (b)  the  measured  steady-state  burning  rate 
as  a  function  of  pressure  and  Initial  temperature,  and  (c)  the  following 
algebraic  relationship. 


♦s  “ 


pr 


(T  -  T  ,) 
ps  pi 


(9) 


which  represents  the  surface  heat  flux  for  steady-state  burning.  Detailed 
justification  for  this  procedure  can  be  found  In  Ref.  2. 
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STEADY-STATE  EXPERIMENTAL  DATA 


1/T 

ps 

Figure  2a 


Figure  2c 

Fig.  2  Schematic  drawing  showing  the  procedure 
to  obtain  the  heat  feedback  function 
from  steady-state  data 
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2.1.2  Numerical  Scheme 


The  governing  partial  differential  equation  for  the  solid 
phase,  Eq.  (3),  was  solved  by  using  an  Implicit  finite  difference  tech¬ 
nique  with  a  variable  mesh  system.  The  variable  grid  spacing  provides 
finer  grid  spacing  near  the  surface.  The  central-difference  method  was 
used  to  approximate  both  time  and  spatial  derivatives.  When  T^  represents 
T  (  Ja  t ,  y^) ,  the  finite  difference  approximation  for  the  time  and  spatial 
derivatives  are: 


ay 


aT 

T J+1  _  tj 
i  i 

(10) 

at 

1 

At 

aT 

j 

bl  Ti-1  +  b2  Ti  +  b3  Ti+1 

(ID 

1 

yl  "  yi-l 

-a2_T- 

2 

j 

a 

c,,  yi  — , 

(12) 


where 


b 


1 


1+f 

2+f 


_ 1 _  (13) 

(1+f)  (2+f) 


c 


1 


S 


2 

2+f  * 


2+f  c  _  1 

1+f  »  J  1+f  (14) 


The  factor,  f,  is  used  to  increase  the  mesh  size  in  the  y  direction: 


yi+ryi 


Ayt  =  (i+f)  ay^j 


(15) 
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Using  the  above  finite  difference  approximations,  the  transient 
heat  conduction  equation  can  be  expressed  in  the  following  form: 


T^1-  TJ 


At 


*2 


Ij _ .(»!  (tJ!}  b3(T^|  +  TJ+l)l 

2  Ay 


ri-l 


a  c 
pr  1 


(16) 


2  )  -  c  (T 

2  (Ayi_i)  KVi  T  i-r  2  t 


+TJ  )  +  =  3<T£!  +  Ti+1)I 


3+14.TJ 

i 


It  should  be  noted  that  the  subscript  pr  for  temperature  has  been  dropped 
here  for  brevity.  Eq.  (16)  was  linearized  by  approximating  the  burning 
rate  at  the  intermediate  time  step,  with  the  average  value  calculated 
between  the  previous  time  value  and  the  current  time  value  in  the  last 
iteration,  1.  e. , 


J+l 


1  ,  r  J 

2  '  bi 


+  r^»k) 

bi 


(17) 


for  the  (k  +  l)th  iteration. 


The  following  expression  is  used  at  the  solid  boundary: 


a  T 

*  y 


j+i 


-  t 


(18) 


where  is  obtained  through  the  burning-rate  map,  details  of  which  are 
given  in  the  next  section.  The  left-hand  Tide  of  Eq.  (18)  is  approximated 
by  a  three-point  backward  difference  .oethod5  and  is  represented  by 


a  t 
Ty' 


j+l 


-p2  (T ,J+1  -  xf1)  +(lW  (T?1  -  T^1) 


_L 


(l  +p) 


where  p  -  l/(l+f ). 
For  the  last  nodal  ooint  (i 


imax). 


_L 


(19) 


imax 


Pi 


(20) 
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The  resulting  set  of  simultaneous  algebraic  equations  was  solved  by  using 
a  standard  tridiagonal  matrix  inversion.  The  propellant  surface  tempera¬ 
ture  thus  obtained  was  used  in  the  Arrhenius  expression,  Eq.  (7),  to 
obtain  the  new  burning  rate.  This  burning  rate  wa;.  then  used  to  obtain 
from  the  burning-rate  map.  The  temperature  profile  was  recomputed, 
based  upon  the  new  value  of  +  ;  before  advancing  to  the  next  time  step, 
the  cycle  was  continued  until  the  solution  converged. 


2.1.3  Implementation  of  the  Dynamic  Burning-Rate  Program 

2. 1,3. 1  Development  of  Burning-Rate  Map  for  Propellant  A 


The  Zeldovlch  burning-rate  map  was  constructed, 
udng  the  steady-state  burning  rate  data  supplied  by  hr.  C.  F.  Price  of 
NWC  (see  Appendix  I).  To  increase  the  number  of  data  points  at  different 
pressures  and  initial  temperatures,  additional  intermediate  data  were 
generated  through  interpolation  to  facilitate  map  reading.  Fig.  3  gives 
the  burning  rate  versus  initial  propellant  temperature  for  different 
pressures.  An  average  temperature  sensitivity  coefficient  of  the  burning 
rate  was  also  obtained  from  the  data.  The  temperature  sensitivity  is 
needed  to  extrapolate  the  available  experimental  data  for  obtaining  the 
propellant  Initial  temperature  below  the  range  of  data  supplied  by  NWC. 
The  average  value  of  the  temperature  sensitivity  coefficient  was  found 
to  be 


9ln  r 
It 


pi 


*  ) 
p 


0.00132  K 


-1 


(21) 


The  burning  rate  was  plotted  against  1/T  for 

ps 

various  pressures.  As  seen  in  Fig.  4,  one  line  fits  all  the  data  well, 
regardless  of  the  pressure.  The  functional  expression  for  this  line 
was  obtained  in  the  Arrhenius  form  given  by 


E 

rb  “  As  6XP  (~  rV  ) 

u  ps 
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Initial  Temperature,  Tp^,  K 

3  Steady-state  burning  rate  for  various  initial 
propellant  temperatures  and  pressures 


NWC  TP  6193 


g 

where  A  ■  2.7126  x  10  cm/sec 
8 

■  34,541  cal/gm-mole 
Ry  «  1.9872  cal/gm-mole“K 

According  to  the  procedure  outlined  In  Fig.  2,  a  plot  of  ♦  (P,  r.  )  was 

obtained,  as  shown  in  Fig.  5.  Unfortunately,  it  is  nearly  impossible  to 

interpolate  *  g  for  different  values  of  r^  and  P  because  of  the  steepness 

of  the  isobarlc  lines  and  che  scarcity  of  the  data  points.  Accordingly, 

this  plot  is  not  used,  Instead,  Fig.  3  and  Eq.  (22)  are  used  together 

with  Eq.  (9)  to  obtain  ♦  g.  An  interpolation  routine,  based  upon  the 

Altken-Lagrange  interpolation  method,  was  implemented  to  Interpolate  data 

(plotted  on  Fig.  3)  for  obtaining  an  equivalent  initial  propellant 

temperature  T  ,  .  Detailed  procedures  for  computing  the  dynamic  burning 

P  eq 

rate  are  given  in  the  next  section. 

2. 1.3. 2  Dynamic  Burning-Rate  Subroutine. 

After  the  onset  of  ignition,  the  burning-rate  sub¬ 
routine  is  used  to  compute  the  Instantaneous  burning  rate.  First,  the 
values  of  the  burning  rate  and  propellant  surface  temperature  are  assumed. 
The  burning-rate  map.  Fig.  3,  is  Interpolated  to  obtain  the  equivalent 
initial  propellant  temperature.  The  heat  feedback,  ♦  ,  Is  calculated 
from  Eq.  (9).  The  transient  heat  conduction  equation  is  solved  by  using 
a  new  value  of  4g  to  obtain  a  new  propellant  temperature  profile  and 
surface  temperature.  T^.  A  new  burning  rate  is  determined  from  the 
Arrhenius  expression,  Eq.  (22).  The  burning-rate  map  is  used  again  to 
obtain  new  equivalent  initial  propellant  temperature.  This  procedure 
is  repeated  until  the  solution  converges.  The  converged  burning  rate  is 
the  desired  dynamic  burning  rate.  A  general  layout  of  the  dynamic 
burning-rate  subroutine  is  given  in  Fig.  6. 
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Fig.  5  Zeldovich  map  constructed  from  steady-state 
burning  rate  data 
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FIG. 


Fluw  Chart  of  the  Computation  Procedure  used  in  Dynamic  Burning 
Rate  Program 
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During  the  calculation,  if  the  burning  rate  goes 
below  the  value  corresponding  to  that  it  274.7  K,  the  equivalent  Initial 
temperature  Is  obtained  by  using  the  temperature  sensitivity  coefficient 
o  . 


P 


(T 

Fi 


) 

eq 


274.7 


In 


{  [rb<P)I274.7  K) 

l  TT^pT] -  / 


(23) 


2. 1.3. 3  Direct  Versus  Indirect  Methods  for  Incorporating  the 
Dynamic  Burning-Rate  Subroutine  Into  CCC 


The  dynamic  burning-rate  routine  and  the  burning-rate 
map  were  coupled  directly  into  the  crack  combustion  code  (CCC).  Several 
runs  were  made  with  this  coupled  program.  The  following  general  observa¬ 
tions  about  the  coupled  program  were  made: 


1.  Depending  upon  the  temperature  profile  at  the 
onset  of  gasification,  extinction  may  immediately 
follow,  due  to  the  shallow  subsurface  temperature 
profile. 

2.  3 P/3  t  has  a  very  strong  influence  on  burning- 
rate  augmentation. 

3.  In  order  to  obtain  a  convergent  solution,  the 
time  increment  for  the  coupled  program  must  be  an 
order  of  magnitude  smaller  than  that  for  CCC  alone. 
This  results  In  a  substantial  increase  of  the 
computer  time  for  each  run. 


In  view  of  the  substantial  increase  In  computation  time,  it  was  concluded 
to  be  more  pragmatic  and  economical  to  develop  a  burning-rate  augmentation 
function.  Results  obtained  from  a  systematic  variation  of  parameters  (such 
as  3 P/3 1,  propellant  subsurface  temperature  profile,  etc.)  can  be  used  to 
obtain  such  a  function. 
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In  the  development  of  burning-rate  augmentation  function,  the  effect 
of  subsurface  temperature  profile  at  the  onset  of  ignition  has  to  be  con¬ 
sidered.  A  reasonable  representation  of  the  temperature  profile  would  be 
an  exponential  form.  Therefore,  the  following  exponential  distribution  of 
the  subsurface  temperature  profile  was  used: 


T 

_££ 

T 

ps 


(24) 


It  should  be  noted  that  when  C  =  1,  the  temperature  profile  reduces  to 
that  of  a  steady-state  burning.  When  C  <  1,  the  temperature  profile  will 
be  more  shallow  than  that  of  the  steady-state  profile.  Therefore,  if 
C  <  1  at  the  onset  of  ignition,  no  dynamic  burning  is  assumed  to  take 
place.  In  general,  at  the  onset  of  ignition,  the  subsurface  temperature 
profile  will  not  be  exponential.  A  least  square  curve-fit  procedure  was 
therefore  employed  to  obtain  the  equivalent  exponential  temperature  pro¬ 
file  at  the  onset  of  ignition.  In  the  development  of  the  burning-rate 
augmentation  function,  C  was  used  as  one  of  the  parameters.  The  burning- 
rate  augmentation  factor  is  defined  as  the  ratio  of  the  dynamic  burning 
rate  at  a  given  pressure  to  the  steady-state  burning  rate  at  the  same 
pressure. 


2.1.4  Calculated  Results 


The  numerical  solution  of  the  dynamic  burning  effect  was  studied 
by  two  methods:  (1)  by  using  a  dummy  main  program  with  a  specified  pres¬ 
surization  rate,  initial  propellant  temperature  profile,  and  initial 
pressures;  (2)  by  directly  coupling  the  dynamic  burning-rate  subroutine 
into  the  crack  combustion  code.  As  mentioned  in  section  2.1.3. 3,  the 
time  increment  for  the  coupled  program  must  be  at  least  one  order  of  magni¬ 
tude  smaller  than  that  for  CCC  alone;  therefore,  for  the  most  part,  the 
first  method  was  used  to  study  the  dynamic  burning  effect. 
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Fig.  7  Is  a  typical  graph  showing  the  effect  of  pressurization 
rate  on  dynamic  burning  rate  augmentation  factor,  r^/r^,  for  a  given  pres¬ 
sure  and  steady-state  initial  temperature  profile.  The  pressurization  rate, 
was  assumed  to  be  constant,  with  pressure  Increasing  raonotonically. 

The  reasons  for  considering  this  type  of  pressurization  process  are:  (1) 
during  the  flame-spreading  process  in  a  solid  propellant  crack,  the 
pressurization  rate  at  a  given  axial  location  remains  more  or  less  con¬ 
stant;  (2)  a  single  variable  (l.e. ,  the  average  value  of  ^  )  can  be 
used  to  study  the  effect  of  pressurization  rate  in  the  parametric  study; 
and  (3)  during  a  DDT  process,  the  pressure  In  the  system  usually  increases 
monotonlcally. 

A  similarity  between  the  shapes  of  various  curves  under  dif¬ 
ferent  pressurization  rates  can  be  seen  in  Fig.  7.  The  burning-rate 
augmentation  factor  starts  from  1  at  t  ■  0.  As  time  increases,  the  augmen¬ 
tation  factor  also  Increases,  reaching  a  plateau.  As  time  increases 
further,  the  curve  passes  through  a  point  of  inflection,  and  finally 
reaches  a  runaway  condition.  The  time  required  to  reach  the  runaway  con¬ 
dition,  or  the  vertical  asymptote,  decreases  as  pressurization  rate  is 
increased.  In  other  words,  the  plateau  (flat)  region  of  the  curve  decreases 
as  the  pressurization  rate  is  increased,  and  finally  disappears  when  the 
pressurization  rate  Is  greater  than  10^  atra/s.  It  is  also  evident  that 
for  pressurization  rates  above  1(T  atra/s,  the  dynamic  burning  rate  is 
important  even  for  short  time  durations.  This  implies  that  the  dynamic 
burning-rate  effect  may  contribute  significantly  to  the  DDT  process,  which 
usually  occurs  at  very  high  pressurization  rates  (»  10^  atm/s). 

Fig.  8  is  similar  to  Fig.  7,  except  that  the  initial  temperature 
profiles  at  the  onset  of  ablation  are  different.  Figure  8  corresponds  to 
the  case,  with  higher  initial  thermal  energy  storage  than  shown  in  Fig. 

7.  By  comparing  these  two  figures,  the  effect  of  initial  temperature 
profiles  can  be  readily  seen.  A  thicker  Initial  temperature  profile  (which 
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Fig.  7  Calculated  dynamic  burning  rate  augmentation 

factor  versus  time  for  C  -  1.0  and  P.  =  82.5  atm 
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corresponds  to  larger  values  of  parameter  C)  results  In  an  increased 
dynamic  burning  rate.  However,  the  time  required  to  reach  the  runaway 
condition  remains  essentially  unaffected  by  the  value  of  C.  Fig.  9  is  a 
plot  of  burning-rate  augmentation  factor  versus  time  for  C  -  1.  5  and 
initial  pressure  of  50  atm.  Comparing  Figs.  8  and  9,  one  can  see  that  as 
the  Initial  pressure  is  decreased,  the  general  nature  of  the  dynamic 
burning-rate  curves  remain  constant;  however,  the  curve  for  a  given  pres¬ 
surization  rate  and  initial  temperature  profile  shifts  to  the  left,  i.e. , 
f >  ••  time  to  reach  the  runaway  condition  decreases.  This  is  in  agreement 
with  the  generally  observed  pressure  effect  on  dynamic  burning.* 

It  should  be  noted  that  the  pressure-time  history  has  a  strong 
influence  on  the  predicted  dynamic  burning  rate.  Numerical  calculations, 
with  a  amp  pressurization  followed  by  a  constant  pressure  period,  show 
that,  the  burning  rate  declines  after  reaching  a  peak,  and  approaches  the 
stea  -state  burning  rate.  These  results  are  also  consistent  with  those 
rep..  d  in  previous  dynamic  burning  studies. 

The  effect  of  the  initial  temperature  profile,  which  is  thinner 
than  the  ,teady-state  profile,  has  also  been  studied.  It  was  found  that 
when  the  initial  temperature  profile  is  shallower  than  the  steady-state 
case,  i.  ...  ,  C  <  1,  the  dynamic  burning-rate  augmentation  during  a  pres¬ 
surization  process  is  much  less  pronounced.  Under  high  pressurization 
rates,  the  burning  rate  may  first  reach  a  peak  followed  by  a  valley  and 

then  increase  again  after  a  long  period  of  time.  At  lower  pressurization 

rates,  the  burning  rate  may  first  reach  a  peak  followed  by  a  monotonic 

decline  until  extinction  is  reached.  Therefore,  the  dynamic  burning  affect 

is  not  taken  into  account  for  cases  with  C  <  1  in  the  development  of  the 
dynamic  burning-rate  augmentation  function. 

2.1.5  Development  of  Burning-Rate  Augmentation  Function 

Several  sets  of  data  were  generated  by  using  the  dynamic 
burning-rate  subroutine  to  facilitate  the  development  of  a  burning-rate 
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augmentation  function.  The  parameters  used  to  generate  the  dynamic 
burning  rate  were:  the  pressurization  rate,  dP/bt;  the  Initial  pressure, 
P^;  and  the  variable  C,  which  governs  the  shape  of  the  Initial  temperature 
profile  [see  Eq.  (24)].  For  the  sake  of  generality,  the  burning-rate 
augmentation  function  was  developed  through  the  use  of  a  set  of  dimension¬ 
less  parameters.  The  reference  parameters  used  for  nondimensionallzation 
are 
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(26) 


where  r°  Is  the  steady-state  burning  rate  obtained  from  Saint-Robert' 8 
type  of  turning  rate  law: 


where  a  ■  0.9359  x  10  ^  (cra/s)/(g/cm2)n 
n  -  0.4108 


(27) 


The  dimensionless  parameters  are 


PR  = 


:£)/(£)* 


(28) 


T  =  t/t* 


(29) 


BR  = 


where  P 

i 


50  atm 


Data  were  generated  for  eight  pressurization  rates,  dp/dt  ■  0,  104, 
10^,  2  x  10^,  5  x  10^,  10^,  and  102  atm/s;  three  initial  propellant 


(30) 

(31) 

5  x  104, 
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temperature  profiles,  C  ■  1,  1.5  and  2;  and  three  Initial  pressures, 

P .  -  34,  50,  and  82.  5  atm.  Since  the  nature  of  the  curves  is  very  similar, 

it  Is  possible  to  obtain  s  suitable  function  which  can  approximate  all  of 

the  curves.  Because  the  runaway  time  or  the  vertical  asymptote  of  the 

burning  rate  for  a  given  PR  and  PI  is  nearly  independent  of  C,  the  fast* 

rising  portion  of  the  curves  in  the  later  stage  can  be  approximated  by 


BR  - 


(32) 


where  a^,  a2,  a3,  and  a^  are  functions  of  PI  only.  In  order  to  obtain  a 
better  fit  for  the  entire  region,  the  following  functional  form  was  used: 


Kt  2 

BR  -  t  1  (^  +  k3*  +  k^r  )  +  a3 


'<  ^  J 

PR 


(33) 


where  the  functional  dependence  of  k's  are 


\  -  k,(PI,  C) 
k2  -  k2(PI,  PR) 
k3  -  k3(PI,  PR) 
kA  -  k4(PI,  PR) 


(34) 


The  coefficients  a^  and  a^  were  approximated  by  second-order  polynomials 
of  PI,  and  are  given  by 


al 

'0.975  -0.6303  1.6235 

'  1 

a2 

0.743  -0.0384  0.1314 

PI 

a3 

“ 

0.465  -0.  1967  0.091 

7 

PI 

a. 

0.5867  0.  1283  0,2781 

4  J 
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The  coefficients  k's  are  given  by 

kj  -  (0. 424  +  0.514  PI  +  2.6159  PI2)  -  (0.508  +  0.204 
+  (0.128  -  0.1379  PI  +  0.0677  PI2)C2 

PI  +  0.1136  PI2)C 

(36) 

k2  -  (1  -  PI)( 1. 673  -  1.81  PR  +  0.272  PR2)  x  1016 

(37) 

k3  -  (1  -  PI)(1, 639  -  1.762  PR  +  0.2648  PR2)  x  1015 

(38) 

k,  -  (1  -  PIX9.226  -  9.98  PR  +  1.5  PR2)  x  1017 

(39) 

The  coefficients  were  obtained  through  the  use  of  an  established  statistical 
analysis  (SAS)  program. 7  Figures  10  to  18  Bhow  the  comparison  of  the  dynamic 
burning  rates  predicted  by  the  burning-rate  augmentation  function  with  those 
obtained  using  the  dynamic  burning-rate  program.  It  can  be  seen  from  these 
plots  that  the  proposed  dynamic  burning-rate  augmentation  function  represented 
by  Eq.  (33)  can  closely  predict  the  burning  rates  obtained  from  the  finite 
difference  solution  for  the  range  of  pressurization  rates,  initial  pressures, 
and  the  initial  propellant  temperature  profiles  considered  in  this  study. 
Figures  16  to  18  show  the  degree  of  scatter  of  the  burning  rates  calculated 
by  the  two  different  methods  (numerical  vs.  augmentation  function)  around 
the  45°  line.  It  is  quite  clear  that  the  agreement  is  good  for  all  three 
initial  pressures  considered.  Therefore,  it  is  believed  that  the  burning- 
rate  augmentation  function  can  adequately  predict  the  dynamic  burning  effect. 

2.1.6  Limitations  in  the  Use  of  Dynamic  Burning-Rate  Augmentation 
Function 


In  spite  of  the  close  agreement  between  the  burning  rates  pre¬ 
dicted  by  the  augmentation  function  and  those  computed  by  the  finite 
difference  method,  care  should  be  taken  in  using  the  augmentation  function 
beyond  the  range  of  parameters  studied.  It  should  also  be  noted  that  the 
application  of  the  augmentation  function  is  limited  because  of  the  follow¬ 
ing  constraints: 
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Fig.  10  Comparison  of  the  dynamic  burning  rate  data  with  the 
prediction  from  burning  rate  augmentation  function 
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Fig. 11  Comparison  of  the  dynamic  burning  rate  data  with  the 
prediction  from  burning  rate  augmentation  function 
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Time,  t,  ms 

Fig. 14  Comparison  of  the  dynamic  burning  rate  data  with  the 
prediction  from  burning  rate  augmentation  function 
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Dynamic  to  Strand  Burning  Rate  Ratio  Predicted 
by  the  Augmentation  Function,  r^Zr” 

Fig.  16  Comparison  of  the  dynamic  burning  rates  calculated 
by  two  different  methods 
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Dynamic  to  Strand  Burning  Rate  Ratio  Predicted 

by  the  Augmentation  Function,  r,  /rp 

h  b 

Fig.  17  Comparison  of  the  dynamic  burning  rates  calculated 
by  two  different  methods 
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Dynamic  to  Strand  Burning  Rate  Ratio  Predicted 

By  the  Augmentation  Function,  r,/r° 

b  b 


Fig.  18  Comparison  of  the  dynamic  burning  rates  calculated 
by  two  different  methods 
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(a)  The  pressurization  rate  was  assumed  to  be  constant  for 
the  entire  period. 

(b)  In  the  generation  of  augmentation  function,  no  pressure 
cutoff  was  considered. 

(c)  The  Initial  propellant  subsurface  temperature  profile  was 
approximated  by  an  exponential  function.  This  may  not 
always  be  an  accurate  representation  of  the  actual 
temperature  profile. 

(d)  The  function  should  not  be  used  for  cases  with  C  <  1. 

(e)  The  Zeldovich  method,  which  was  used  to  calculate  the 
dynamic  burning  rate,  breaks  down  under  erosive  burning 
conditions.  Therefore,  the  proposed  function  should  not 
be  used  for  situations  involving  erosive  burning. 
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2.2  Technological  Gaps  In  the  Existing  Granular  Bed  Combustion  Models 
for  DDT  Studies 

8~  12 

In  recent  years, four  one-diraensional  two-phase  combustion  models 
have  been  developed  for  prediction  of  the  combustion  processes  in  granular 
propellant  beds.  These  models  were  independently  developed  at  The  Pennsylvania 
State  Univeristy,  Gough  Associates,  Inc. ,  The  University  of  Illinois,  and 
Calspan  Corp.  for  interior  ballistic  predictions.  They  were  reviewed  and 
summarized  in  a  JANNAF  Workshop^  conducted  three  years  ago. 

Some  modifications  and  improvements  in  these  models  have  been  made 
since  the  JANNAF  Workshop.  At  the  Pennsylvania  State  University,  improve¬ 
ments  have  been  made  in  the  following  three  areas*  The  intragranular  stress 

and  particle-wall  friction  measurements  in  granular  bed  were  made  and 

13 

correlations  were  obtained.  The  speed  of  sound  transmitted  through  a 

packed  granular  bed  as  a  function  of  porosity  was  also  obtained.  To 

determine  the  actual  flow  rates  from  percussion  primers,  a  primer  character- 

14 

lzation  study  was  conducted.  The  results  of  this  study  were  used  as  input 
data  for  granular  bed  combustion  calculations.  The  boundary  condition 
treatment  at  the  ends  of  the  granular  bed  has  been  further  improved.  The 
most  up-to-date  procedures  for  boundary  condition  calculations  are  reported 
in  Ref.  15.  The  comparison  of  theoretical  predictions  and  experimental  data 
are  given  in  Refs.  15  and  16. 


The  flow  resistance  measurements  through  granular  propellant  beds 

made  of  cylindrical  grains  were  conducted  by  Robbins  and  Gough ^  and  their 

latest  experimental  results  are  in  close  agreement  with  Kuo-Nydegger's 
18 

correlation.  For  large  caliber  systems,  the  multi-dimensional  effect 

is  important,  and  Gough  has  recently  developed  a  two-dimensional  convective 

19 

flame  spreading  model.  Special  emphases  were  placed  in  the  description 
of  the  internal  boundaries  so  that  the  effect  of  bag  material,  liners  and 
additives  can  be  studied. 
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Krler  has  changed  his  earlier  continuum  mechanics  model  Into 

20 

a  separated  two-phase  flow  model  for  studying  DDT  phenomena  in  granular 
propellant  beds.  A  number  of  empirical  correlations  were  altered  fot 
using  them  beyond  the  ranges  for  which  these  correlations  were  developed. 
The  lack  of  suitable  constitutive  relationships  were  noted. 

Besides  the  constitutive  relationships  for  flow  resistance,  heat 
transfer,  intragranular  stress,  etc.,  there  are  a  number  of  technological 
gaps  that  should  be  narrowed  in  order  to  extend  the  existing  granular  bed 
combustion  models  for  DDT  studies.  These  are  given  in  the  following. 

1.  Better  ignition  criterion  has  to  be  established.  The 
gasification  of  reactive  species  may  not  result  In  an 
instantaneous  reaction  to  generate  heat  at  the  location  of 
gasification.  Also,  the  products  In  the  form  of  condensed 
phases  may  have  non-negllglble  effect  on  ignition.  The 
effect  of  shock  wave  propagation  In  granular  material  can 
also  have  Important  contribution  to  Ignition.  The  shear, 
compression,  and  friction  effects  should  also  be  considered. 

2.  The  evaluation  of  intragranular  stresses  transmitted 
through  aggregates  of  solid  propellants  under  dynamic 
conditions  must  be  performed. 

3.  The  fracture  of  grains  and  resulting  increase  In  the 
specific  surface  area  of  the  propellant  under  dynamic 
compression  should  be  studied. 

4.  Dynamic  and  erosive  burning  effects  during  the  transition 
process  from  deflagratidn  to  detonation  should  also  be 
evaluated. 
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5.  The  treatment  of  extremely  steep  pressure  gradients  near 
the  wave  front  in  the  granular  bed  should  be  Incorporated 
in  the  solution  scheme. 

6.  Understanding  of  the  generation  of  hot  spots  in  solid 
propellants  must  be  advanced. 

7.  Mechanical  properties  and  fracture  behavior  of  the 
propellant  under  extremely  rapid  pressurization  conditions 
must  be  characterized. 
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2. 3  Documentation  of  the  Fixed  and  Mobile  Granular  Bed  Combustion  Codes 
(MGBC  end  FGBC ) 


To  facilitate  the  use  and  understanding  of  the  granular  bed  combustion 

program,  two  user’s  manuals  were  prepared:  one  was  designed  for  mobile 

21  22 
granular  bed  combustion  and  the  other  for  fixed  granular  bed  combustion. 

These  manuals  outline  theoretical  models,  numerical  solution  methods,  and 
use  of  computer  codes  for  predicting  the  physical  processes  involved  in 
mobile  and  fixed  granular  bed  combustion  of  solid  propellants,  respectively. 
The  manuals  provide  details  of  the  following:  (1)  description  of  the 
physical  processes;  (2)  basic  assumptions  used  in  the  model;  (3)  formula¬ 
tion  of  the  model,  using  the  basic  conservation  equations  for  single  or 
two-phase  flow;  (4)  initial  and  boundary  conditions;  (5)  use  of  method  of 
characteristics  for  extraneous  boundary  conditions;  (6)  procedures  for  the 
propellant  surface  temperature  calculations;  (7)  empirical  correlations 
used  in  the  model;  (8)  numerical  scheme  for  the  interior  points;  (9) 
numerical  scheme  for  the  boundary  conditions;  (10)  description  of  the 
main  program  and  the  subroutines;  (11)  input  data  for  the  Initial  and 
restart  runs;  (12)  typical  input  and  output  of  the  program;  (13)  finite- 
difference  equations  for  the  governing  equations;  (14)  procedures  for 
determining  the  flow  properties  at  the  boundaries  of  the  granular  bed; 
and  (15)  glossary  of  important  variables  used  in  the  program.  Both  the 
MGBC  manual  (78  pages)  and  the  FGBC  manual  (40  pages)  have  been  delivered 
to  NWC  for  use  by  NWC  personnel.  It  is  believed  the  manuals  give  suffi¬ 
cient  details  to  enable  users  to  easily  modify  the  program  for  particular 
applications. 
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2. 4  Delivery  of  MCBC  and  FGBC  Programs  to  NWC  and  Consultation  for  Their 

Implementation 

Both  MGBC  and  FGBC  programs  were  delivered  to  NWC  durtng  the  early 
phase  of  the  contract  period.  Before  delivery,  the  programs  were  revised 
by  deleting  unnecessary  portions  and  by  adding  a  large  number  of  comment 
cards  to  facilitate  understanding  of  the  program  by  new  users. 

MGBC  was  delivered  In  the  form  of  a  magnetic  tape,  and  FGBC  ms 
delivered  in  the  form  of  cards.  Sample  case  input  data  and  output  list¬ 
ings  were  also  prepared  and  delivered  with  each  source  program.  An 
updated  version  of  the  MGBC  program  was  mailed  to  Mr.  C.  F.  Price  In 
September  1979.  Two  sets  of  sample  data  with  the  description  of  Input 
Format  were  prepared  for  the  updated  MGBC. 

During  the  contract  period,  consultations  were  provided  for  Imple¬ 
mentation  of  these  programs  onto  the  computer  at  NWC.  These  programs 
have  been  successfully  executed  after  some  modifications  by  MR.  C.  F. 
Price  at  MWC. 
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III.  SUMMARY  AND  CONCLUSIONS 


1.  Dynamic  burning  effects  in  the  combustion  of  solid  propellant 
cracks  were  studied,  using  the  Zeldovich  quasi-steady  flame 
approach. 

2.  Transient  heat  conduction  equation  for  the  solid  propellant  was 
solved  with  the  given  Zeldovich  burning  rate  map  by  using  an 
implicit  finite  difference  numerical  scheme. 

3.  Computed  results  indicate  that  the  dynamic  burning  effect  becomes 
stronger  for  higher  pressurization  rates,  higher  initial  energy 
storage  in  the  solid  propellant  (i.e. ,  thicker  temperature 
profile  at  the  onset  of  ablation),  and  lower  initial  pressure. 

4.  For  monotonically  increasing  pressure  with  a  constant  pressuriza¬ 
tion  rate,  a  runaway  condition  for  the  dynamic  burning  rate  was 
observed.  The  time  required  to  reach  the  runaway  condition 
decreases  as  the  pressurization  rate  is  increased,  or  when  the 
initial  pressure  is  decreased;  the  runaway  time  is  found  to  be 
very  weakly  dependent  upon  the  initial  propellant  temperature 
profiles  considered  in  this  study. 

5.  The  dynamic  burning  effect  was  found  to  be  important  in  the  com¬ 
bustion  processes  in  solid  propellant  cracks.  The  dynamic  burning 
consideration  was  incorporated  into  the  crack  combustion  code  (CCC) 
in  two  ways.  The  first  method,  which  involves  a  direct  coupling 

of  the  dynamic  burning  subroutine  with  CCC,  requires  extremely 
small  time  increments  and  was,  therefore,  not  used.  The  second 
method  involves  the  development  of  a  dynamic  burning-rate 
augmentation  function  for  an  indirect  coupling  with  CCC. 

6.  A  dynamic  burning-rate  augmentation  function  was  generated. 

This  function  fits  the  burning  rate  data,  calculated  from  the 
Zeldovich  method,  over  a  broad  range  of  conditions.  The  augmen¬ 
tation  function  was  incorporated  into  CCC. 
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A  literature  survey  was  conducted  on  recent  studies  In  granular 
propellant  bed  combustion  In  order  to  examine  their  applies  illty 
for  DDT  studies.  Technological  gaps  In  the  existing  granular  bed 
combustion  models  for  DDT  study  were  listed. 

Two  user's  manuals  were  prepared  and  delivered  to  NWC.  One  was 
designed  for  mobile  granular  bed  combustion  (MGBC)  and  the  other 
for  fixed  granular  bed  combustion  (FGBC). 


The  updated  computer  codes  for  MGBC  and  FGBC  were  also  delivered 


NWC  TP  6193 


IV.  REFERENCES 


1.  Kuo,  K.K. ,  and  Coates,  G.  R. ,  "Review  of  Dynamic  Burning  of  Solid 
Propellants  In  Gun  and  Rocket  Propulsion  Systems,"  the  Sixteenth 
Symposium  (International)  on  Combustion,  MIT,  Cambridge,  Massachusetts, 
pp.  1177-1192,  August,  1976. 

2.  Summer f ield,  M. ,  Cavney,  L.H. ,  Battista,  R.A. ,  Kubota,  N. ,  and  Isoda,  H. , 
"Theory  of  Dynamic  Extinguishment  of  Solid  Propellants  with  Special 
Reference  to  Nonsteady  Feedback  Law,”  Journal  of  Spacecraft  and  Rockets. 
Vol.  8,  No.  3,  pp.  551-558,  March,  1971. 

3.  Kuo,  K.K. ,  Chen,  A. T. ,  and  Davis,  T. R. ,  "Convective  Burning  In  Solid- 
Propellant  Cracks,”  AIAA  Journal.  Vol.  16,  No.  6,  pp.  600-607,  June,  1978. 

4.  Kuo,  K.K. ,  McClure,  D. R. ,  Chen,  A.T. ,  and  Lucas,  F.G. ,  "Transient 
Combustion  In  Solid  Propellant  Cracks,"  Naval  Weapons  Center  Report  No. 

NWC  TP  5943,  October,  1977. 

5.  Kuo,  K.K. ,  Kovalcln,  R.L. ,  and  AUcman,  S.J.,  "Convective  Burning 

In  Isolated  Solid-Propellant  Cracks,"  Annual  Report  to  Naval  Weapons 
Center,  February,  1978. 

6.  Kuo,  K.K. ,  "A  Summary  of  the  JANNAF  Workshop  on  Theoretical  Modelling 
and  Experimental  Measurements  of  the  Combustion  and  Fluid  Flow  Processes 
In  Gun  Propellant  Charges,"  13th  JANNAF  Combustion  Meeting,  CPIA 
Publication  281,  Vol.  1,  pp.  213-233,  September,  1976. 

7.  Barr,  A. J.,  Goodnight,  J.H. ,  Sail,  J.P. ,  and  Helwlg,  J.T. ,  "Statistical 
Analysis  Systems, "SAS  Institute,  Inc.,  1976. 


50 


NWC  TP  6193 


8.  Kuo,  K.K. ,  Koo,  J.H. ,  Davis,  T.  R. ,  and  Coates,  G. R» ,  "Transient 
Combustion  in  Mobile  Gas-Permeable  Propellants,”  Acta  Astronautlca, 

Vol.  3.,  No. 7-8,  pp.  574-591,  July,  1976. 

9.  Gough,  P. S. ,  "Fundamental  Investigation  of  the  Interior  Ballistics 

of  Guns,"  Final  Report,  Contract  N00174-73-C-0501,  SCR-R-75,  July,  ] 974. 

10.  Gough,  P. S. ,  "The  Flow  of  a  Compressible  Gas  Through  an  Aggregate  of 
Mobile  Reacting  Particles,"  Ph. D.  Thesis,  McGill  University,  December, 
1974. 

11.  Krier,  H. ,  and  Ra jan,  S. ,  "Flame  Spreading  and  Combustion  in  Packed 
Beds  of  Propellant  Grains,"  AIAA  Paper  75-240,  AIAA  13th  Aerospace 
Sciences  Meeting,  January,  1975. 

12.  Fisher,  E. B. ,  and  Trippe,  A.  P. ,  "Mathematical  Model  of  Center  Core 
Ignition  in  the  175  mm  Gun,"  Calspan  Report  No.  VQ-5163-D-2,  March,  1974. 

13.  Kuo,  K.K. ,  Moore,  B. B* ,  and  Yang,  V.,  "Measurements  and  Correlation 
of  Intragranular  Stress  and  Particle-Wall  Friction  in  Granular 
Propellant  Beds,"  16th  JANNAF  Combustion  Meeting,  Naval  Postgraduate 
School,  Monterey,  California,  September,  1979. 

14.  Kuo,  K.K. ,  Moore,  B. B. ,  and  Chen,  D.Y.,  "Characterization  of  Mass 

Flow  Rates  for  Various  Percussion  Primers,"  7th  International  Colloquium 
on  Gasdynamlcs  of  Explosives  and  Reactive  Systems,  Gottingen,  West 
Germany,  August,  1979. 

15.  Chen,  D.Y.,  and  Kuo,  K.K. ,  “Boundary  Condition  Treatment  of 
Transient  Two-Phase  Flows  in  Gun  Systems,"  AIAA  Paper  No.  80-0080, 
to  be  presented  at  Pasadena,  California,  January,  1980. 


51 


NWC  TP  6193 


16.  Davis,  T.  R. ,  and  Kuo,  K.K. ,  "Experimental  Study  of  the  Combustion 
Process  In  Granular  Propellant  Beds,"  Journal  of  Spacecraft  and  Rockets. 
Vol.  16,  No.  4,  pp.  203-209,  July-August,  1979. 

17.  Robbins,  F. ,  and  Gough,  P.  S. ,  "Influence  of  Length  and  Diameter 

of  Cylinders  on  Packed  Bed  Flow  Resistance,"  16th  JANNAF  Combustion 
Meeting,  Naval  Postgraduate  School,  Monterey,  California,  September,  1979. 

18.  Kuo,  K.K. ,  and  Nydegger,  C.C. ,  "Flow  Resistance  Measurement  and 
Correlation  in  a  Packed  Bed  of  WC  870  Ball  Propellants,"  Journal  of 
Ballistics,  Vol. 2,  No.  1,  pp.  1-26,  1978. 

19.  Gough,  P.  S.  ,  "Two  Dimensional  Convective  Flamespreading  in  Packed  Beds 
of  Granular  Propellant,"  Contract  Report  ARBRL-CR-00404,  July,  1979. 

20.  Krier,  H. ,  and  Kezerle,  J.A. ,  "A  Separated  Two-Phase  Flow  Analysis  to 
Study  Deflagration-to-Detonation  Transition  (DDT)  in  Granulated 
Propellant,"  Seventeenth  Symposium  (International)  on  Combustion. 
pp.  23-34,  1978. 

21.  Kuo,  K.K. ,  and  Kumar,  M. ,  "Program  Description  Manual  for  Mobile 
Granular  Bed  Combustion  (MGBC)  Code,”  Submitted  to  NWC,  China  Lake, 
California,  January,  1979,  (78  pages). 

22.  Kuo,  K.K. ,  Kuo,  O.J.,  and  Kumar,  M. ,  "Program  Description  Manua]  for 
Fixed  Granular  Bed  Combustion  (FGBC)  Code,"  Submitted  to  NWC,  China  Lake, 
California,  March,  1979,  (40  pages). 


52 


NWC  TP  6193 


APPENDIX  I 


STEADY  STATE  BURNING  DATA  FOR  PROPELLANT  A 


Burning  Rate  (cm/s)  of  Propellant  A  as  a  function 
of  pressure  and  initial  temperature 


Initial 
Temp. ,  K 


321.9 


15,470 

36,910 

mm 

0.711 

m 

0.729 

0.516 

0.759 

Pressure,  gf/cm 


72,770 


107,900 

143,400 

1.143 

1.229 

1.179 

1.252 

1.234 

1.  306 

Surface  temperature  at  Ignition  (K)  of  Propellant  A 
as  a  function  of  pressure  and  initial  temperature 


Initial 
Temp. ,  K 


321.9 


Pressure,  gf/cm* 


15,470 

36,910 

72,770 

107,900 

143,400 

773.278 

788.155 

799.150 

805.678 

808.423 

774.809 

789.069 

800.030 

306. 850 

809.128 

776,761 

790.5295 

801.734 

808. 578 

810.738 
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INPUT  DATA  DESCRIPTION  FOR  DYNAMIC  BURNING  RATE  PROGRAM 


Computer  Symbol  Symbol  In  Analysis  Description 


Units 


Card  1  (5D14.7) 

SAENG  E 

a 


ARRHC  Ag 

WTF 

PDYMAX 


PRATE  dP 

dt 

Card  2  (5D12.5) 

DELTAT  t 

ALPHAP  a 

P 

BREXP  n 


A 


a 


Surface  activation  energy 
used  in : 

r,  -  A  exp(-E  /Ru/T  ) 
b  s  a  ps 

Arrhenius  coefficient 


cal/g-mol 


cm/s 


Weight  fraction  in  calcula¬ 
ting  new  assumed  r^.0<WTF-l 

Maximum  pressure  above  which 
no  dynamic  burning  effect  is 
considered 

The  ramp  pressurization  rate 


g/cm 


2 


g/cm  -s 


Time  increment  s 

2 

Thermal  diffusivity  cmVs 

Burning  rate  exponent  in 
Saint-Robert's  Law 

2 

Burning  rate  coefficient  (cm/s)/(g/cm  ) 

in  Saint-Robert's  Law 


YTPDIP 


(yT  ) 

P 


dip 


The  depth  below  t!\e  turning  cm 

propellant  surface  at  which 
the  dipstick  ignition  criterion 
is  applied. 


Card  3  (2D12.5.  IX.  815) 

TPI  T 

Pi 

FAC  FAC 


KX 


NMAX 


Propellant  initial  temperature  K 


Factor  used  in  the  expansion  of 
the  grid  size  normal  to  the 
burning  surface 
Ay (I)  =  (1  +  FAC)Ay(I-l) 

The  integer  to  designate  the 
axial  location  along  the  crack. 
KX»1  represents  the  crack  entrance 
mesh  point;  KX*NMAX+1  represents 
the  closed  end  of  the  crack. 

The  total  number  of  Ax  intervals 
along  the  crack. 
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Computer  Symbol  Symbol  in  Analysis 
Card  3  (2D12.5.  IX,  815)  (con’t) 

ITRMAX 

NDTMAX 

IRSTRT 


NPRINT 


I0VRWT 


NEXTIN 


Card  4  (5D12.5) 
CTS1 


CTS2 


TREFL 


T 


ref 


Description 


The  maximum  number  of  itera¬ 
tions  allowed  for  the  conver¬ 
gence  of  burning  rate. 

The  maximum  number  of  At  in 
an  initial  or  restart  run. 

The  flag  to  indicate  the  re¬ 
start  run. 

IRSTRT=1,  for  restart  run 
IRSTRT=0,  for  initial  run 

The  flat  to  control  the  print¬ 
ing  of  output.  Print  after 
every  NPRINT*At . 

A  flat  to  control  the  over¬ 
write  option  of  the  subsurface 
temperature  profile.  I0VRWT^O, 
an  exponential  temperature  pro¬ 
file  is  used  to  overwrite  the 
initial  temperature  profile. 

A  flag  to  control  the  burning 
rate  to  avoid  extinction. 
NEXTIN^O,  instantaneous  burning 
rate  cannot  be  less  than  the 
steady-state  value. 


Coefficient  1  used  in  the  temp¬ 
erature  sensitivity  expression: 


a  =  C 
p  Ts 


+  C„  P 
1  Ts2 


Coefficient  2  used  in  the  temp¬ 
erature  sensitivity  expression: 

°p  =  CTS1  +  CTS2P 


The  lowest  available  T^  data 

used  as  a  reference  *  temper 
ature  for  the  equivalent  Tp 
through  the  use  of  temper-  * 
ature  sensitivity. 
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Units 


1/K 


cm2/g-K 


K 
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